d 

o 


RESEARCH  AND  DEVELOPMENT  BRANCH 
DEPARTMENT  OF  NATIONAL  DEFENCE 

CANADA 


CO 


ESEARCH 


USHMENT  OTTA\ 

i  i  %  y 

DRLO  TECHNICAL  NOTE  81-18 


A  COMPUTER  PROGRAM  TO  GENERATE 


CL. 

O 

C  "S 


Lul 
— J 

ul 


RAIViAN  SPECTRA 


JANUARY  \m 
OTTAWA 

ea  05  13  os® 


PROJECT  HO.  256 


RESEARCH  AND  DEVELOPMENT  BRANCH 


DEPARTMENT  OF  NATIONAL  DEFENCE 
CANADA 


DEFENCE  RESEARCH  ESTABLISHMENT  OTTAWA 

DREO  TECHNICAL  NOTE  81-18 


A  COMPUTER  PROGRAM  TO  GENERATE 
THEORETICAL  COHERENT  ANTI-STOKES 
RAMAN  SPECTRA 


by 

T.  Parameswaran 
Roy  Ball  Associates  Ltd. 
and 

D.R.  Snelling 

Energy  Systems  Section 
Energy  Conversion  Division 


!  fy.orovd  tor  poMto 

rwrtributtoB  UaMwMe* 

r  i  — ' - 


JANUARY  1982 
OTTAWA 


UNCLASSIFIED 


ABSTRACT 


A  computer  code  for  generating  and  plotting  theoretical  CARS  (Coherent 
Anti-Stokes  Raman  Spectroscopy)  spectra  is  described.  This  code,  which  was 
obtained  from  Dr.  A.C.  Eckbreth  of  the  United  Technologies  Research  Centre, 
was  modified  for  use  on  a  Sigma-9  computer.  The  theory  and  computational 
procedure  involved  in  generating  CARS  spectra  is  discussed  and  the  input 
parameters  and  variables  used  in  the  code  are  described. 


Theoretical  nitrogen  spectra,  which  are  presented  for  a  range  of 
experimental  conditions  are  discussed. 
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Description  d'un  code  machine  permettant  de  creer  et  de  representer 
graphiquement  des  spectres  theoriques  obtenus  par  la  spectroscopie  Raman 
"Anti-Stokes"  coherente.  Ce  code  a  ete  fourni  par  le  Dr.  A.C.  Eckbreth  du 
United  Technologies  Research  Centre  et  a  ete  modifie  pour  permettre  son 
emploi  avec  l'ordinateur  Sigma-9.  II  y  a  aussi  la  description  de  la  theorie 
et  de  la  procedure  de  calcul  qui  entrent  dans  la  creation  de  ces  spectres, 
ainsi  que  des  paramStres  et  des  variables  d' entree  employes  dans  ce  code. 

On  y  trouve  enfin  la  representation  des  spectres  theoriques  de  l'azote 
presentes  pour  toute  une  gamoe  de  conditions  experimentales. 
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INTRODUCTION 


A  project  to  develop  Coherent  Anti-Stokes  Raman  Spectroscopy  (CARS) 
as  a  combustion  diagnostic  tool  was  recently  initiated  at  DREO.  This  technique 
promises  to  provide  the  capability  of  non-perturbing  measurements  of  temperature, 
pressure,  and  composition  (major  species)  in  real  combustion  environments.  CARS 
can  provide  good  spatial  (<1  mm^)  and  temporal  (<10  ns)  resolution  in  high 
interference  environments  (such  as  sooting  flames  which  are  highly  luminescent). 

In  order  to  determine  the  temperature  and  composition,  experimental 
spectra  are  fitted  to  those  derived  theoretically.  Therefore,  concurrently 
with  the  development  of  the  CARS  apparatus,  we  are  undertaking  development  of 
the  theory  and  software  necessary  for  the  generation  of  CARS  spectra.  This 
work,  in  part,  is  being  undertaken  by  Roy  Ball  Associates  Ltd.  under 
contract  2-ST81-000055. 

We  have  been  greatly  aided  by  having  a  computer  code,  developed  by 
Dr.  A.C.  Eckbreth's  group  at  United  Technologies  Research  Centre  (UTRC), 
made  available  to  us. 

This  code  was  developed  by  R.J.  Hall  at  UTRC^>^)  on  the  UNIVAC 
1110  Computer.  We  have  suitably  modified  the  above  program  so  that  it  can  be 
run  on  the  Sigma-9  computer  at  CRC/DREO.  The  plots  of  the  resulting  spectra 
are  obtained  by  the  additional  use  of  PLOT-10  software  with  the  4015-1 
Tektronix  Display  terminal.  In  this  note  we  will  explain  the  theoretical 
principles  underlying  this  computer  program  and  then  discuss  some  of  the 
procedural  details  of  the  Fortran  code. 


THEORY 


The  theory  of  coherent  anti-stokes  Raman  spectroscopy  has  been  dealt 
with  extensively  (3»7).  In  brief,  CARS  is  a  nonlinear  optical  phenomenon 
and  involves  three  wave  mixing.  When  two  laser  beams  of  frequencies  u)£  and 
u interact  in  a  medium,  three  wave  mixing  produces  a  resultant  coherent  beam 
with  the  frequency  2u)£  -  and  this  is  the  CARS  beam.  The  mixing  occurs 
for  all  samples  but  the  intensity  of  the  CARS  signal  is  greatly  enhanced  when 
u)£  -  ws  approaches  a  Raman  frequency  of  the  medium. 

The  interaction  of  the  laser  radiation  with  the  medium  occurs  through 
the  third  order  nonlinear  electric  susceptibility  denoted  by  x^.  gives 

rise  to  an  induced  polarization  field  which  acts  as  a  source  term  in  Maxwell's 
wave  equation.  On  solving  the  wave  equation  one  gets  the  expression  for  the 
power  of  the  CARS  signal,  as 
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Thus  the  CARS  power  P3  varies  linearly  with  the  Stokes  power  Ps,  quadratically 
with  the  pump  power  P^  and  is  also  proportional  to  the  square  of  the  third 
order  susceptibility.  has  resonant  and  non-resonant  contributions. 

x(3)  _  yR  +  (2) 


X«)  can  be  obtained  semiclassically  or  quantum  mechanically  and  is  given 
by 


R  2Nc4 

do 

x  ‘ 

d  Si 

j 

_  id2  -(«d4  -  ids)2  -i  Ti  (id£-ws)  _ 

where  N  is  the  number  density,  (dj  is  the  Raman  frequency,  Ij  is  the  Raman  line 
width  and  do/dft  is  the  scattering  cross  section.  It  is  clear  that  becomes 
large  when 


=  U>£-Ws 


(4) 


If  we define  the  frequency  detuning  as 


A  id  =  -  (u)£-(ds)  and  if  Aid/idj  <<  1 

1  s  1 


OJj-Aw  Wj 


1  +  Aus 


*3  j 


(5) 


(6) 


Then  eqn  (3)  becomes 


r  2Nc* 
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r 

da 

i 

dQ 

j 

2  Au  ir 

{  i-  i  \ 

(7) 


t  -  v  2Nc4  A,  da 

Let  Kj  m*S3T  j  dfl 


rj _1 


(8) 


where  Aj  is  the  population  difference  between  the  upper  and  lower  states 
of  the  rotation-vibration  transitions  with  frequency  id j .  Summing  over 
all  such  wj  (with  AJ  =  0,±2)  that  lie  near  u>£  -  ids 
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X  (3> 


l 
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2AViri 


(9) 


The  Raman  cross  section 


do 

dQ 


do 

3TZ 


T 


can  be  expressed  as  follows  for  the  Q  branch. 


2Mu>o 


2 

a 


4s 


(10) 


where  M  is  the  reduced  mass,  (jJq  the  angular  frequency  of  the  molecular 
oscillator,  b*l  are  the  Placzek-Teller  coefficients,  a  is  the  derivative  of 
the  mean  molecular  polarizability  with  respect  to  the  inter  nuclear 
coordinate  and  y  the  similar  derivative  of  anisotropy  and  v  is  the  vibrational 
quantum  number  of  the  initial  level. 

X(3)  can  thus  be  calculated  from  equation  (9)  Then,  treating  N  and 
Xm  as  parameters  the  CARS  power  can  be  obtained  at  a  given  temperature  T. 

In  the  above  analysis  the  pump  frequency  ii)£  and  the  Stokes  frequency 
ojg  have  been  assumed  to  be  ideally  monochromatic.  In  practice  the  pump 
laser  may  have  a  finite  width  and  a  broadband  Stokes  laser  can  be  used  to  get 
multiplex  CARS  spectra  in  a  single  shot.  Assuming  the  modal  phases  to  be 
uncorrelated,  the  laser  power  can  be  written  as 


P«, 


(<*>£>  dW£ 


Ps  = 


Ps  (ws )  d  s 


(11) 


and  the  CARS  power  (u>)  can  be  obtained  by  convoluting  over  the  laser 
widths. 


f  2 

P3  {<*>}  «\»  P£{a)£}dW£  P£(u)-6}  Ps{oo£-6}  j  x{6  }  |  d6 


(12) 


Generally,  Aw^  <<Aws  an<*  very  small.  Therefore  Pgtwj-g}  can  be 


replaced  by  Ps  {u)£i  J  -6}  and  the  integral  can  be  simplified  to 
P3(w}  'V  P^  Ps  {w-6}  |  x!<5}|2d6 


(13) 


The  computer  program  developed  by  R.  Hall  utilizes  the  above  convolution 
integral  to  compute  the  CARS  intensity  as  a  function  of  frequency.  A  Gaussian 
distribution  is  assumed  for  the  spectral  densities  P^w^),  Ps(u)g) 


i.e. 


P£  {U)£} 

Pn  {  to  } 
s  1  S 


exp 


exp  ^ 


width 


Au) 
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width 
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'  2 
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(14) 

(15) 
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where  ■  u >-  u)£,Au>s  *  w-u)s  and  width  denotes  the  laser  line  width. 

The  CARS  power  obtained  from  equation  (12)  should  be  again  convolved 
over  the  detector  slit  function,  which  is  usually  taken  to  be  triangular. 
The  final  CARS  power  thus  found  can  be  compared  with  experimental  data. 


GENERATION  OF  THEORETICAL  CARS  SPECTRA 
ENERGY  LEVELS  AND  RAMAN  FREQUENCIES 

For  a  given  temperature  T,  the  main  program  calculates  the  energy  levels 
and  rotational  -  vibrational  frequencies  of  a  diatomic  molecule,  based  on 
the  Boltzmann  distribution.  For  a  vibrating  rotator  the  term  values  are 
given  by  (H) 

T'  =  =  G  (v)  +  Fv  (J)  (16) 

where  G(v)  is  the  vibrational  and  FV(J)  the  rotational  contributions. 


G(v)  =  we(v+i)  -  uje  xe (v+i ) 2 

+  weye  (v+i)3  +  t oeze  (v+£)4 


(17) 


where  v  is  the  vibrational  quantum  number  and  u>e,wexe,u>eye ,  weze  are 
spectroscopic  constants. 


FV(J)  *  ByJ  (J+l)  -  DyJ2  (J+l)2  (18) 

By,  Dv  are  found  from  the  Dunham  expansion. 

By  *  Be  -  oie(v+£)  +ye  (v+i ) 2  (19) 

Dy  =  De  +  3e  (v+i)  (20) 

where  again  ae,  3e,  Ye  are  spectroscopic  constants. 

Knowing  the  energy  levels  E(v,J)  the  transition  frequencies  Wy.j  can 
be  canclulated. 

The  number  of  molecules  in  a  vibration  level  of  energy  Ey  is  Ny  * 
e*p(-Ey/KT)  -  exp  (-1.439  G(v)/KT). 
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Relative  to  the  ground  state  v=0  the  term  values  are 

G0  (v)  =  u)e(v)  -  u)exe  (v2+v)  + -  (21) 

Thus  when  exp  (-1.439  GQ(v)/KT)  =  exp  (-1.439  wev/KT)  (22) 


is  small  the  number  of  electrons  in  the  level  is  small.  When  this  quantity 
is  $0.01  the  corresonding  v  +  1  is  taken  as  the  limit  of  vibrational  levels 
that  contribute  significantly  to  the  CARS  intensity. 

The  number  of  molecules  in  the  rotational  level  J  of  a  rigid  rotator 
is  given  by 


Nj  «  (2J+1)  e"  BVJ  (J+l)  /KT 


(23) 


Since  (2J+1)  increases  linearly  with  J,  Nj  goes  through  a  maximum  before 
going  to  zero.  _ 


Setting  9Nj/3J  =  0  we  get  Jmax  = 


KT 


2Bvhc 


-i 


for  B  *  10.44  cm  1,  T  =  300°  K,  Nj  -»■  0  around  3.5  times  (see  Fig. 

Ref.  (11). 


9 


Therefore 

the  upper  limit  . 
contribute  to  x  . 


KT 

2Bvhc 


is  used  in  the  program  as  a  criterion  for  finding 
of  the  number  of  rotational  states,  that 


POPULATION  FACTORS 

The  vibrational  partition  function  is  determined  in  Subroutine  QVIB. 


Qvib  =  l  exP  ("G(v)/KT)  (24) 

v 

This  subroutine  also  computes  the  vibrational  population  factors 
exp(-G(v)iKT)  =  /N 
Qvib 

where  Ny  is  the  number  of  molecules  in  state  v  and  N  is  the  total  number 
of  molecules. 

Similarly  function  QROT  computes  the  rotational  partition  function 

Qrot  =  l  <2J+1>  exP  <-F(J)/KT)  (25) 

J 
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Subroutine  BOLTZ  then  calculates  the  combined  Boltzmann  population  factor 
required  in  the  expression  for  the  susceptibility  •  The  population 
factor  Ay,j  is  given  by 


v*J 


-hcFL(J)  v 
e  KT 


-he  FU(J ) 
e  KT 

Qrot 


NvU 

N 


(26) 


where  L  denotes  the  lower  state  and  U  the  upper  state. 


NUCLEAR  STATISTICAL  WEIGHT  FACTOR 

Molecular  wave  functions  are  symmetric  or  antisymmetric  under 
exchange  of  nuclei.  Generally  even  J  functions  are  symmetric  and  odd  J 
antisymmetric  or  vice  versa.  When  the  nuclear  spin  1=0  the  symmetric 
and  antisymmetric  states  never  combine  and  we  see  missing  lines  in  the  Raman 
rotational  spectrum  (e.g.  0£).  When  nuclear  spin  /  0  these  spatially 
symmetric  and  antisymmetric  levels  combine  but  with  different  statistical 
weights. 

Consider  1^  =  I2  =  i 
Then  total  I  =  1,  0 

I  =  1+3  spin  functions  that  are  symmetric 
I  =  0*1  spin  function  that  is  antisymmetric 

Therefore  even  I  spin  functions  pair  with  spatially  symmetric  functions  and 
odd  I  spin  functions  pair  with  spatially  antisymmetric  functions. 

Thus  the  statistical  weight  factor  is  3  for  odd  J  and  1  for  even  J 
if  even  J  states  are  symmetric  and  odd  J  states  are  antisymmetric,  and 
vice  versa.  As  a  second  example, 

Consider  1^  =  I2  =  1 
Then  I  =  2,  1,  0 

I  =  2~*5  symmetric  spin  functions 
I  =  l-*3  antisymmetric  spin  functions 
I  =  0*1  symmetric  spin  function 

Again  even  I  states  pair  with  even  J  (if  symmetric)  and  odd  I  states  pair 
with  odd  J  (if  antisymmetric)  and  vice  versa.  For  nitrogen  even  J  states  are 
symmetric  and  odd  J  states  are  antisymmetric;  therefore  the  nuclear  spin 
degeneracy  factor  is 

6  for  even  J 
3  for  odd  J 

These  nuclear  statistical  weight  factors  are  used  in  summing  over  the  various 
rotational  levels. 
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RAMAN  SCATTERING  CROSS  SECTION 

In  its  present  form  this  program  computes  the  susceptibility  for  Q 
branches  only.  For  Q  branches  the  scattering  cross  section  is  given  by 


do 


1  r2a! 

=  l—s- 

L.  c 

4 

2Mw 

'Q  l 

o 

{v+1} 


described  in  eqn  (9).  Knowing  the  factors  in  question,  do/dft  can  be  easily 
found.  The  Placzek-Teller  factors  are  from  ref  (10) 


J  _  J  (J+l) _ 

J  (2J-1)  (2J+3) 


NON-RESONANT  SUSCEPTIBILITY 


The  non-resonant  part  of  y  is  real  and  proportional  to  the  number 
density  , 


i.e. 


=  Net' 


where  a'  is  the  molecular  polarizability.  N  the  number  density  at  temperature 
T  and  pressure  P  is  given  by 

P  x  y*1*  at  NTP 

The  number  density  N  of  a  gas  in  a  mixture  of  gases  is  given  by 


Xmole  x  pressure  x  1.013  *  10^ 
1.38  x  10-16  x  T 


!7) 


where  \iole  is  the  molar  fraction  of  the  gas  in  question.  The  factor  1.013  x 
106  comes  from  changing  atmosphere  to  dynes/cm^. 


PRESSURE  BROADENING  OF  RAMAN  LINES 

The  CARS  code  calculates  the  Raman  line  width  according  to  the 
equation 


r  (J)  =  P. (8T-0,71-1E.6T  ”1,45j) 


(28) 


where  T  is  the  temperature,  P  the  pressure,  J  the  rotational  quantum  number 
and  r(J)  is  the  line  width  FWHM  in  cm”^  for  900K  ^T$  2400K.  This  value  of 
r ( J )  is  based  on  the  isolated  line  approximation. 
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The  parameters  in  the  above  equation  were  obtained  by  Hall  by  a 
least  square  fit  of  N2-Q  branch  line  widths  obtained  by  Owyoung  and  Rahnd^) 
with  the  theoretical  line  width  model  predicted  by  Bonamy  and  Robert(^). 

CORRECTION  FOR  x(3^ 

This  program  computes  the  anharmonicity  correction  and  the  correction 
due  to  centrifugal  distortion  that  contribute  to 

These  are  obtained  from  the  vibrational  and  vibration-rotation  matrix 
elements  of  the  dipole  moment  operator  with  anharmonic  terms  in  the  potential 
(15,16). 


The  anharmonicity  correction  is  obtained  from 

B-i  {v'>v+1^  (v+1)^  (4-6 (v-H)  d  +  ll(v+l)>  b2 

£  (4-6d  +  lib2) 


(29) 


where  the  quantities  d,b,  will  be  described  in  the  next  section. 

The  correction  due  to  centrifugal  distortion  is  given  by  the  matrix 
element,  below  for  Q  branches. 


v,J-»v+l,J 

X1 


where  s 


-  (v+3)r2+i (v+l)r2+  £(13v  +  16)rb 

-  £  (13v+10)rb 

=  3J  (J+D^i  ;  r  =  J  (J+l)y3/2 
2 


(30) 


NUMBER ICAL  SCHEME  AND  CONVOLUTION  INTEGRALS 

Thus  far  we  have  discussed  only  the  frequency  independent  terms  of 
The  frequency  appears  in  the  denominator  of  x^  in  the  detuning  factor. 

2Au)j  =  2  {coj  —  oos  >  } 

The  numerical  scheme  in  the  main  program  selects  a  frequency  range  A,  divides 
it  into  a  suitable  number  of  discrete  intervals,  and  generates  the  frequency 
shifts  (W£-us)  «X^  is  then  calculated  at  all  these  points  summing  over  all 
the  states  in  the  Q  branch. 

The  numerical  scheme  also  divides  the  chosen  frequency  range  A  into 
intervals  at  which  the  CARS  power  P3  ^  is  calculated.  P3  (w)  is  calculated 
by  numerical  evaluation  of  the  convolution  integral  in  eqn.  (13). 
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Simpsons  rule  for  numerical  integration  according  to  which 

|  ydx  =  yo  +  4^1  +  2Y2  +  4?3  + -  +  y2in  (31) 

*o 


for  an  even  number  of  intervals,  is  used  to  evaluate  the  integral. 

Subroutine  SLIT  uses  the  same  integration  rule  for  convolution  over 
a  triangular  slit  function.  The  effect  of  varying  the  slit  width  can  be 
studied  by  calling  SLIT  several  times. 

Plots  of  the  CARS  spectra  computed  by  the  program  are  obtained  with 
the  Tektronix  terminal  by  using  PLOT-10  software. 

R.  Hall  has  included  additional  codes  to  compute  corrections  to  the 
susceptibility  arising  from  collisional  line  narrowing  (17).  This  effect 

becomes  significant  at  high  pressures  i.e.  >  20  atmospheres.  We  have  not  used 
this  part  of  the  code  as  of  now.  The  details  of  this  section  will  be  described 
in  a  later  note. 


DESCRIPTION  OF  INPUT  PARAMETERS  AND  VARIABLES  IN  THE  COMPUTER  CODE 

1.  WE,  XWE,  YWE,  ZWE 

These  are  the  molecular  spectroscopic  constants. 

2.  BE,  DE,  Bl,  B2,  D1 

are  also  spectroscopic  constants  with  the  following  relation  with 
Herzbe'"g's  0-1)  notation. 

Bl  ae 
B2  -  Ye 
D1  -*  ge 

3.  SNUC1,  SNUC2  are  the  statistical  weights  for  even  and  odd  levels 
respectively,  due  to  the  nuclear  spin  degeneracy. 

4.  CONST  is  the  constant  term  in  given  by  ^j_0  where  M  is  the  reduced 

mass  and  wo  the  angular  frequency  of  the  molecular  oscillator. 

5.  ALPHA,  GAMMA  are  the  derivatives  with  respect  to  the  internuclear  distance 
of  the  mean  molecular  polarizability  and  anisotropy,  respectively. 

6.  AWID,  AWIDE,  BWID,  BWIDE 

are  parameters  used  to  find  the  temperature  and  rotational  quantum 
number  dependent  pressure  broadened  line  width  WIDTHP(J)  (Refer  to 
R.  Hall's  paper  in  App.  Spcy(^)). 

7.  T  is  the  mean  rotational  -  vibrational  temperature  -  an  input  parameter. 

8.  PRESS  is  the  gas  pressure. 
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9.  XMOLE  is  the  molar  fraction  of  the  species. 

10.  CHINR  Is  the  nonresonant  part  of  x  and  is  also  an  Input  parameter. 

11.  0MEGA1,  WID1,  OMEGA2,  WID2 

are  the  pump  laser  frequency,  width,  stokes  laser  frequency  and  width 
In  that  order. 

12.  NS,  WS1IT 

are  the  number  of  slit  functions  and  the  slit  widths  for  those  functions. 

13.  W3MIN,  W3MAX 

define  the  limits  of  the  spectral  range  for  calculations  of  • 

14.  NV  is  the  number  of  vibrational  states 
NJ  is  the  number  of  rotational  states. 

15.  DENSY  is  the  number  density. 

16.  W3PTS+1  is  the  number  of  points  at  which  the  convoluted  CARS  power  is 
calculated. 

17.  Al,  A2,  GAMM,  BA,  DA 

are  needed  for  finding  the  anharmonicity  correction  VCOR  and  the 
centrifugal  distortion  VJCOR  to  X^  (15,  16).  In  relation  to  Bouanich 
and  Brodbeck's  article 

BA  -+b;  DA  -*d;  Aj-*-a^;  GAM  -*y 

y  =  2BE/WE  b  =  Al  J  GAMM/ 2 

d  =  A2  GAMM/ 2 

18.  PTQ  (J)  are  the  Placzek-Teller  factors. 

19.  WIDTHP  (J)  is  the  temperature  and  rotational  number  dependent  pressure 
broadened  Raman  line  width. 

20.  E(I,J)  are  the  energy  levels 

LAMQ  (I,J)  are  the  transition  frequencies  for  the  Q  branch. 

21.  WTEST  ■  LAMQ  (I  J)  -  DOMEGA  is  the  frequency  detuning  factor. 

22.  WPTS+1  is  the  number  of  frequency  points  for  X^  calculation. 

23.  FREQ  (K+l)  stores  the  above  Raman  frequencies. 

24.  0MEGA3  (K+l)  stores  the  frequencies  defined  by  W3PTS  for 
convolution  integration. 

25.  POWER  (K+l)  stores  the  CARS  power  at  points  defined  by  Freq  (K+l). 
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26.  P0WER3  (K+l)  is  the  array  of  CARS  power  values  after  convolution  over 
laser  line  width. 

27.  P0WER4  (K+l)  is  the  similar  array  after  slit  convolution. 

28.  INITT,  BIN ITT,  CHECK,  D SPLAY,  XFRM,  YFRM,  XTYPE,  YTYPE,  TINPUT, 

FINITT  are  subroutines  from  PLOT-10  software. 


In  the  original  version  of  the  program  these  plot  routines  are  included 
in  the  main  section.  We  have  chosen  to  separate  the  plotting  package.  We 
calculate  and  store  the  required  arrays  and  plot  them  later  as  required,  with 
the  plotting  program  CARSPLOT  consisting  of  the  above  PLOT-10  routines. 


RESULTS  AND  DISCUSSION 


The  computer  program  described  earlier  was  used  to  generate  the 
theoretical  CARS  spectra  of  nitrogen.  These  spectra  are  displayed  in  Figures 
1  through  3. 

These  spectra  were  calculated  with  a  pump  laser  FWHM  of  0.8  cm-!  and 
a  width  of  175  cm"l  for  the  Stokes  laser.  The  values  of  non-resonant 
susceptibility,  the  molecular  polarization  parameters  and  the  spectroscopic 
constants  of  nitrogen  used  in  the  input  data  were  the  same  as  those  used 
by  R.  Hall  (1). 

Figure  1  depicts  the  N2  CARS  spectrum  at  atmospheric  pressure  and 
at  a  temperature  of  2100°K.  In  Fig.  la  the  Q-branch  spectrum  of  N2  for 
monochromatic  incident  beams,  is  shown.  The  stronger  band  is  the  0-1 
fundamental  and  the  weaker  one  corresponds  to  the  1-2  hot  band.  The 
spectrum  is  composed  of  two  distinct  envelopes  with  different  intensities. 

For  nitrogen,  the  nuclear  spin  is  1  and  this  leads  to  a  statistical  weight 
factor  of  6  for  the  even  J  rotational  levels  and  a  factor  of  3  for  the  odd  J 
levels.  Therefore  the  even  J  levels  have  an  intensity  about  four  times  that 
of  the  odd  levels  and  correspond  respectively  to  the  strong  and  weak  spectral 
envelopes. 

In  Fig.  lb  the  CARS  power  after  convolution  over  the  finite  widths  of 
the  laser  beams,  is  shown.  Gaussian  spectral  profiles  have  been  assumed  for 
both  the  pump  and  the  Stokes  lasers. 

The  resulting  CARS  power  should  be  convolved  once  more  over  the 
instrumental  slit  functions  before  comparing  with  an  experimental  CARS 
signature.  The  convolved  spectrum  for  a  slit  function  of  1.0  cm"l  is 
shown  in  Fig.  lc. 

The  present  program  can  perform  this  slit  convolution  (using  a  triangular 
slit  function)  for  any  chosen  number  of  slit  functions.  In  Figs.  2a  through 
2c  we  see  the  convolved  spectra  for  the  slit  functions  0.8,  1.0  and  2.7  cm-1 
respectively.  It  is  clear  that  0.8  and  1.0  cm"*  yield  high  resolution  spectra 
while  2.7  cni"!  leads  to  lower  resolution.  The  prominent  peak  noticed  in  the 
hot  band  of  the  higher  resolution  spectra  is  due  to  a  spectral  overlap  between 
a  high  J,  0  -*■  1  branch  and  lower  J,  1  -*■  2  branch. 
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The  effect  of  temperature  on  the  theoretical  CARS  spectra  is  shown  in 
Fig.  3.  A  low  resolution  slit  width  of  2.7  cm~^  was  chosen  for  this  purpose. 
Temperatures  used  were  from  1200°  -  2400°K  in  steps  of  300°K.  The  hot  band 
is  sensitive  to  the  temperature,  and  since  the  rotational  structure  is  not 
resolved,  it  can  be  used  for  temperature  measurement.  The  0  ■+•  1  band 
however  is  not  greatly  affected  by  temperature. 
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FIGURE  CAPTIONS 

Fig.  1.  Q-branch  theoretical  CARS  spectrum  of  ^  at  T  =  2100°K, 

P  *  1  atm. 

Fig.  la.  CARS  spectrum  for  monochromatic  frequencies. 

Fig.  lb.  The  spectrum  convoluted  over  finite  laser  widths. 

Fig.  lc.  The  spectrum  after  slit  function  convolution,  slit  width  -  1.0  cm~l 

Fig.  2.  Effect  of  slit  functions;  a)  0.8  cm~l 

b)  1.0  cm-l 

c)  2.7  cm"l 


Fig.  3.  Variation  of  N2  CARS  spectrum  with  temperature  slit  width 
=  2.7  ctn-1,  P  =  1  atm. 
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